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Deep 610-MHz Giant Metrewave Radio Telescope observations of 
the Spitzer extragalactic First Look Survey field - III. The radio 
properties of Infrared-Faint Radio Sources 
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ABSTRACT 

Infrared-Faint Radio Sources (IFRSs) are a class of source which are bright at radio frequen- 
cies, but do not appear in deep infrared images. We report the detection of 14 IFRSs within 
the Spitzer extragalactic First Look Survey field, eight of which are detected near to the lim- 
iting magnitude of a deep 7?-band image of the region, at 7? ~ 24.5. Sensitive Spitzer Space 
Telescope images are stacked in order to place upper limits on their mid-infrared flux den- 
sities, and using recent 610-MHz and 1.4-GHz observations we find that they have spectral 
indices which vary between a = 0.05 and 1.38, where we define a such that Sv = SqV^", and 
should not be thought of as a single source population. We place constraints on the luminos- 
ity and linear size of these sources, and through comparison with well-studied local objects 
in the 3CRR catalogue demonstrate that they can be modelled as being compact (< 20 kpc) 
Fanaroff-Riley Class II (FRII) radio galaxies located at high redshift (z ~ 4). 
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1 INTRODUCTION 

Infrared-Faint Radio Sources (IFRSs) are a class of radio-bright or 
infrared-faint objects id entified in the Aus tralia Telescope Large 

', Area Survey (ATLAS: iNorris et al.J |200^ . The ATLAS survey 
consists of deep 1.4-GHz radio observations over two of the 
southern Spitzer Wid e-area Infrared Extragalactic Survey (SWIRE; 

, iLonsdale et alj2003i) fields, which have deep infrared data between 
3.6 and 16 /xm publicly available. The Spitzer Space Telescope 
dWerner e t al. 2004) infrared observations of the SWIRE fields are 
sensitive enough that it was expected that all radio sources de- 
tected in the ATLAS survey which were in the local U niverse would 
be visible in the SWIRE images ( iNorris et al.ll2006l) . whether the 
radio emission was due to star-formation or an Active Galactic 
Nuclei (AGN). However, it was found that some radio sources 
did not have an infrared counterpart at any of the Spitzer wave- 
lengths. A total of 55 IFRSs of varying flux densities h ave been de- 
tected within the Chandra Deep Field South (CDFS: lNorris et afl 



2006) and European Large- Area ISO Survey South- 1 (ELAIS-Sl; 
Middelberg et al. 2008) fields, with several having 1.4-GHz flux 
densities of greater than 5 miy, and the brig htest having a 1 .4-GHz 
flux density of 26.1 mIy dNorris et all200^ . 

The possibility that these sources are obscured star-forming 
galaxies with infrared flux densities falling just below the sensitiv- 
ity limit of Spitzer would seem to be ruled out by the production of 
stacked infrared images centred on the radio source positions of 22 
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IFRSs dNorris et al.l200^ which show no evidence for any infrared 
counterparts. Norris et aU ( |2007|) made observations of two of their 
bright IFRSs using Very Long Baseline Interferometry (VLBI), and 
found that one of them had a compact core with angular size less 
than 0.03 arcsec (a linear size below 260 pc at any reasonable 
redshift), although they were unable to detect the second target. 
They concluded that the radio emission is driven by AGN activity 
rather than star-formation, making IFRSs either high redshift radio- 
loud quasars, or ab normally highly obs cured radio galaxies at more 
moderate redshifts. INorris et al.l ( |2007|) considered the spectral en- 
ergy distribution of the radio-loud quasar 3C273, and showed that 
placing it at low redshift could not reproduce both the detection at 
1.4 GHz and non-detection at 3.6 ptm for the source with the com- 
pact core, but moving 3C273 to z = 7 would fulfil the requirements. 
However, the lack of information at other wavelengths makes it dif- 
ficult to determine exactly what type of source the IFRSs are. These 
sources are selected through being radio-bright, making the use of 
radio diagnostics important for determining their type. 

The first of the large-scale surveys carried out by the Spitzer 
Space Telescope was the Spitzer extragalactic First Look Survey 
(xFLS). Four square degrees, centred on 17*'18™00', +59°30'00" 
(12000 coordinates, which are used throughout this work), were ob- 
served with two instruments - the Infrared Array Camera (IRAC; 
Fazio et al. 2004) at 3.6, 4.5, 5.8 and 8 Mm, and the Mu ltiband 
Imaging Photometer for Spitzer (MIPS; iRieke et al]|2004h at 24, 
70 and 160 p,m. This data is comparable in depth to the SWIRE 
survey, with catalogue completeness limits (5cT, apart from the 3.6- 
jlm band which has its completeness limit set at 7(7) of 20, 25, 
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100 and 100 /xJy in the four IRAC bands. There are two deep ra- 
dio sur veys of the region, at 1.4 GHz with the Very Large Array 
(VLA; ICondon et aLll2003h and at 610 M Hz with the Giant Me- 
trewave Radio Telescope (GMRT; iGarn et al. 2007, hereafter Pa- 
per I). These have comparable resolution (5 arcsec^ at 1.4 GHz; 
5.8 X 4.7 arcsec-^ at 610 MHz) and sensitivity (23 /iJy beam^' at 
1.4 GHz; 30 /xJy beam^^ at 610 MHz), and cover approximately 
the same region of sky as the infrared obser vations. An R-hw A sur- 
vey of the xFLS field has been carried out jpadda et al.ll2004h with 
a limiting Vega magnitude of R = 25.5, and estimated 50 per cent 
completeness limit of R = 24.5. The deep and complementary in- 
frared, optical and radio observations make this region a good area 
to search for other IFRSs, and the availability of multi-frequency 
radio data will allow the spectral index of any IFRS to be ob- 
tained, providing a much greater understanding of the source prop- 
erties. 

In section |2] we describe the selection criteria for forming a 
sample of Infrared-Faint Radio Sources. We consider the avail- 
able radio and infrared properties of these sources in section [3] 
create stacked infrared images of the sources in order to test 
whether the sources lie slightly below the Spitzer detection thresh- 
old, and place upper limits on their infrared flux between 3.6 and 
160 /Jm. We calculate the radio spectral index of each source, 
and demonstrate that IFRSs are made up of flat, steep and ultra- 
steep spectrum sources. In section |4] we place limits on the lumi- 
nosity and linear size of these sources, and compare well-studied 
sources fr om the Revised Revised Third Cambridge (3CRR) cat- 
alogue of iLaing. Rilev & Longaid jl983l) to the observed IFRS 
population. We demonstrate that it is possible to model the flat, 
steep and ultra-steep-spectrum IFRSs with the spectral energy 
distribution (SED) of a local radio galaxy, placed at high red- 
shift, and conclude that it is possible to explain the IFRS popu- 
lation without requiring new, exotic source types. A flat cosmol- 
ogy with the best-fitting parameters from the Wilkinson Microwave 
Anisotropy Probe ( WMAP) five-year data of Q./^ = 0.74 and Hq = 
72 km s^' Mpc^^ dPunklev et al.l2008l) is assumed throughout this 
work. 



2 SAMPLE SELECTION 
2.1 Initial radio source sample 

The xFLS region was observed at 1 .4 GHz bv lCondon et~ZH2003l) 
with the Very Large Array (VLA), covering ~ 4 deg'^ with a noise 
level of 23 flJy beam^'. We created a sample of radio sources at 
1.4 GHz within the region covered by the VLA a nd GMRT surveys, 
using Source Extractor (SExtractor; iBertin & Arnoutj 
Il99d) and the technique described in Paper I. A peak flux den- 
sity of greater than 0.5 mjy was required for all sources, equiva- 
lent to a signal-to-noise ratio of ~ 20. This ensures that the sources 
are unquestionably real, and that sources which are found to be 
infrared-faint will be extreme examples of their class - the inclu- 
sion of fainter radio sources in the sample would reduce the signifi- 
cance of the lack of an infrared detection. Within the region covered 
by the VLA and GMRT surveys, 511 radio sources were detected. 
We reject 107 of these for being located outside the coverage area 
of one or more of the infrared images, leaving 404 1.4-GHz radio 
sources. 



where we define the radio spectral index a such that Sy = SqV 




(a) A double radio source, with (b) A radio source with slightly 
infrared counterpart lying extended emission covering an 

between the two radio lobes. infrared counterpart. 



Figure 1. Two examples of radio sources with catalogued infrared counter- 
parts which were not initially classified as being a match. Contours are from 
the 610-MHz data, and are 0.5, 1, 2, 4, 8 mJy beam"', while the grey-scale 
shows the 3.6-/im image, and ranges between ± the 3.6-flm noise level of 
2.9 ^Jy pixer'. 

A similar catalogue of radio sources was created from the 610- 
MHz image, requiring sources to be detected at the 5c7 level but 
with no further selection criteria. At the edges of the 610-MHz 
image, 5cT corresponds to a flux density of --^ 750 ^ujy beam"', 
decreasing to ~ 150 /iJy beam"' towards the pointing centres. 
Sources at the 1.4-GHz peak flux limit will therefore be seen any- 
where in the 610-MHz survey if they have a spectral index a 
greater than 0.5, so some flat-spectrum 1.4-GHz sources may not 
be assigned 610-MHz counterparts using this criterion. 

The majority of the 1.4-GHz sources (362/404; 90 per cent) 
had a 610-MHz counterpart within 1.5 arcsec, equivalent to one 
pixel in the radio images. A further 24 had counterparts within 
3 arcsec. The remaining 18 sources were inspected in order to iden- 
tify why no 610-MHz counterpart was found. Seven were either one 
component of a double radio source, or extended radio sources with 
610-MHz counterparts > 3 arcsec away; two were close to a bright 
source and were concealed by the increased noise surrounding it 
(see Paper I for details on the effects of bright sources on the 610- 
MHz image); the remaining nine were near to the edge of the 610- 
MHz mosaic, and not detected due to the increased noise due to the 
decreasing gain of the GMRT primary beam. All 1 8 sources had an 
infrared counterpart (see section 12. 2t , and were rejected from the 
sample, leaving 386 sources to make up the initial 'radio-bright' 
sample. 

2.2 Removing infrared counterparts 

We follow 'Norris et alj ilOO ^ in rejecting sources from the radio- 
bright sample if they have an infrared counterpa rt within 3 arcsec 
in any of the publ i shed Spitzer source catalogues l lLacv et al.l200^ ; 
iFadda et alJl2006l ; lFraver et alj|2006) . This generated an initial list 
of 94 IFRS candidates, which required visual inspection in order 
to check for the presence of infrared counterparts which had not 
been identified as being related to the radio sources. Composite 
infrared-radio maps were generated for each source, after rotating 
the geometry of the IRAC images to match the radio observations. 
The majority of the candidate sources (49/94; 52 per cent) were def- 
inite components of double or triple radio galaxies with an infrared 
counterpart present, but not co-located with the radio emission - 
Fig-Et shows an example of one of these sources, with 610-MHz 
contours overlaid on the 3.6-flm image. The counterpart is suffi- 
ciently far from the centres of the radio components that it was 
not classified as a match, but is clearly associated with the dou- 
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(a) A 'clear' uncatalogued (b) A 'possible' uncatalogued 

counterpart. counterpart. 



Figure 2. Two radio sources with uncatalogued infrared counterparts. Con- 
tours are from the 610-MHz data, and are 0.5, 1, 2, 4, 8 mJy beam^', while 
the grey-scale shows the 3.6-/Jm image, and ranges between ± the 3.6-/im 
noise level of 2.9 /i Jy pixel" ' . 

ble radio galaxy. Four sources were located close to very bright 
infrared sources, which would have concealed any true infrared 
counterparts. There were 21 radio sources which should have been 
matched to a source in the IRAC catalogue, due to either having ex- 
tended radio emission which covered an infrared source, or having 
a slightly greater separation between the infrared and radio source 
centres than 3 arcsec due to the lower resolution of the radio im- 
ages - Fig.[Tj5 shows an example of one of these sources. One ra- 
dio source showed no evidence for infrared emission, but was in a 
noisy region of the 610-MHz image, and the measurement of the 
radio flux was very uncertain. After rejecting all of these types of 
source, there were 19 remaining radio sources with no counterparts 
within the infrared catalogues. 

Five of these sources showed clear evidence for an uncata- 
logued infrared counterpart in the 3.6-/xm image, with four of these 
also showing evidence for an uncatalogued counterpart in the 4.5- 
/xm image. Fig. |2ji shows an example of one of the radio sources 
with a clear 3.6-;tim counterpart, which was not included in the 
IRAC catalogue. We reject these five sources from the sample. A 
further six sources show possible evidence for an uncatalogued in- 
frared counterpart, with a slight increase in infrared emission being 
seen at or close to the centre of the radio emission. We retain these 
sources in the sample, but identify them as having possible counter- 
parts in the remainder of this analysis - we will test for a statistical 
increase in infrared flux at the centre of each radio source in sec- 
tion [3TT] There are 14 IFRSs remaining in the sample, which are 
shown in Fig. |3] All images are 60 x 60 arcsec^ in order to look 
for extended structure. The grey-scale ranges between ±C7 for the 
3.6-/im image noise level of CT = 2.9 /iJy pixep' in order to be 
sensitive to faint infrared emission. The 1.4-GHz contours, and 4.5, 
5.8 and 8-^m grey-scale have a similar morphology. The properties 
of the IFRSs are given in Table[T] and throughout this work, desig- 
nations of the IFRSs will be taken from the 610-MHz FLSGMRT 
catalogue from Paper I. 

2.3 Comparison with previous samples 

Two p revious samples of IFRSs have been constructed. lNorris et al.l 
1 I2OO6I) found 22 IFRSs in the Australia Telescope (AT) observa- 
tions of the 3.7 deg^ CDFS field, while Middellserg et alj bOOSl) 
found 3 1 IFRSs in the observations of the 3.9 deg^ ELAIS-S 1 field. 
The source density of bright IFRSs (with I.4-GHz flux greater than 
I mJy) is comparable in both studies to the sample presented in 
this work, which has a density of ~ 2.5 deg~^, although since we 
have set a lower limit to the radio flux for our sample selection. 



we find fewer faint sources. While large-area and relatively shal- 
low radio surveys such as the I.4-GHz Faint Images of th e Radio 
Sky at Twenty-cm (FIRST jBecker. White & Helfandll 19951) survey 
could be used to find bright IFRSs, the requirement for deep in- 
frared observations means that future studies of this type of source 
are still likely to be limited to smaller well-studied regions such as 
the xFLS and SWIRE fields. 

The number of radio sources which are detected 
increases rapidly with decreasing flux density (e.g. 
ISevmour. McHardv & Gurii] |2004|) . with the radio sky being 
dominated by AGN sources above ~ 1 miy at 1.4 GHz, and the 
contribution from star-forming galaxies becoming important below 
--^ 1 mJy . The fact that we do not see a rapid increase in the 
number of IFRSs at lower radio fluxes is a further indication that 
these sources are unlikely to be obscured star-forming galaxies. 

The re solution of the AT radio images, at 11x5 ar csec^ 
jNorris et a l. 2006) or 10 x 7 arcsec^ (Middelb erg et alj [2008^, is 
poorer than either the VLA or GMRT observations used in this 
work (5 arcsec^ and 5.8 x 4.7 arcsec^ respectively), while the 
SWIRE infrared observations used in the CDFS and ELAIS-S 1 
studies are slightly more sensitive than the xFLS observations. 
While both of these effects may slightly alter the selection biases 
behind the different IFRS samples, w e have followed a si milar se- 
lection procedure to that described in lNorris et al. I j2006l) and be- 
lieve that the sources identified in the three works should be com- 
parable. 



3 INFRARED-FAINT RADIO SOURCE PROPERTIES 
3.1 Source stacking 

A 'stacked' infrared image will reveal whether the IFRS population 
is made up of sources that have infrared fluxes which are slightly 
below the detection threshold. If this is the case, the stacked images 
will show a faint infrared source co-located with the radio source 
centres. A stacking experiment implicitly assumes that the sources 
being stacked are all of the same type - in section [331 we demon- 
strate that this is not the case, with the IFRS population being made 
up of sources with flat, steep and ultra-steep radio spectra. How- 
ever, stacked images will still give an indication of whether there is 
some infrared flux present which is just below the detection limit 
for i ndividual sources. 

iNorris et alj (|2006^ looked at the mean flux in a series of 
stacked IRAC images, and found no detections of IFRSs. How- 
ever, it has been shown (e .g. lWhite et alj200'7l ; lBeswick et al.l20o3 : 
iGam & Alexandeill200 8h that a median stacked image is more reli- 
able than a mean image, since the median image is more robust to 
the presence of a few outlier pixels. We created 14 'cut-out' infrared 
images with size 60 x 60 arcsec^, centred on the radio source posi- 
tions of the IFRS show n in Fig. [3] The in frared background varies 
slightly across the field ('Lacv et alj20o3) and was calculated using 
SEXTRACTOR and subtracted from each cut-out image. Each pixel 
of the stacked image was then created by calculating the median 
value of that pixel from the 14 cut-out images. 

Figs. |4^ to d show the stacked images for the four IRAC 
bands, with the grey-scale varying between ±ct/\/T4 for the rel- 
evant IRAC noise level CT (where tr = 2.9, 5, 20, 20 /iJy for the 
four bands). There is an increase in flux density seen in the centre 
of the 3.6- and 4.5-/lm images, while the 5.8- and 8-/im images 
show no evidence for emission. This is confirmation that at least 
some of the potential uncatalogued infrared counterparts are real. 
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(a) J171141. 3+591813 (b) J171141. 8+591910 (c) J171141.9+591726 (d) J17 1149.8+591727 (e) J171315.5+590302 




(f) J171404.0+595317 (g) J171436.6+594456 (h) J171453.8+594329 (i) J171634.2+601258 G) J171909.8+585346 




(k) J171928.4+583927 (1) J172046.5+584729 (m) J172057.6+601558 (n) J172315.5+590919 



Figure 3. 3.6-^m grey-scale images of the 14 IFRSs, with 610-MHz radio contours overlaid. The designation of each source from the 610-MHz FLSGMRT 
catalogue is given below each subfigure. The grey-scale ranges between ±f7, where a is the noise level of 2.9 /iJy pixel"' , in order to accentuate faint infrared 
emission. The radio contours are at 0.5, 1, 2, 4 and 8 mjy beam"' . The resolution of the radio image is shown by the ellipse in the lower left of each image. 



Table 1. The Infrared- Faint Radio Sources within the Spitzer extragalactic First Look Survey field. Column 1 gives the designation of each source from the 
610-MHz catalogue ofi Gam et al ]j2007), column s 2 and 3 give the 1 .4-GHz and 610-MHz flux densities, column 4 gives the radio spectral inde x g, colunm 5 
gives the R-band magnitude of each source from lFadda et alj i2004h . column 6 gives the deconvolved angular diameters of each source from lCondon et all 
)2003i) and column 7 identifies those sources which have been classified as having potential uncatalogued infrared counterparts. 
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Figs. to h show the stacked images created by using only the 
eight IFRSs which were not identified as having potential uncata- 
logued infrared counterparts - the lower number of sources used 
for the stacking makes these images look noisier, but there is still a 
slight indication of an increase in infrared flux density in the centre 
of the 3.6- and 4.5-/xm images. 

In order to account for the possibility of a detection near to the 



noise level, we place 3ct/\/8 upper limits on the infrared flux den- 
sities of these sources of 3. 1 flly at 3.6 fxra, 5.3 flly at 4.5 flm, and 
21.2 /ijy at 5.8 and 8 /xm, and set a 24-/im limit of 210 /iJy (equal 
to the 24-flm catalogue flux limit for the main survey region), and 
limits of 10 and 50 mJy for the 70-/xm and 160-/im data, again 
from the catalogue flux density limits. 

Star-forming galaxies are observed to follow a tight correla- 
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(e) Ji.6-p.m image, N = S. (f) 4.5-/Jm image. A' = 8. (g) 5.8-/im image, N = S. (h) 8-/im image, N = 8. 

Figure 4. Stacked IRAC images. The grey-scales for each image correspond to the range ±a/\/N, for comparison with Fig. |3] where a = 2.9, 5, 20 and 
20 /I Jy pixel^ ' for the four bands. The radio source locations are at the centre of each stacked image. Images (a) to (d) show the stacked images in the four 
IRAC bands from all 14 IFRSs, while (e) to (h) use only the 8 IFRS from Table[T]which have not been classified as having a possible uncatalogued infrared 
counterpart. 



tion between their mid- or far-infrared flux density, 5ir, and their 
1.4-GHz flux density, 5i .4, which can be quantified by the logarith- 
mic flux density ratio q^^ = lo g]r\(Sm /S] a) (e.g. lAppleton et aU 
l2004lGam. Ford & Alexandeil boOS. hereafter Paper II). The non- 
detections of IFRSs in the 24 and 70-flm Spitzer bands allows us 
to place upper limits of ^24 < —0.7, and ^70 < 1 - these values are 
significantly below the ty pical values seen for s tar-forming galaxies 
of (724 ~ 1 ™d qjQ ~ 2 jAppleton et al.ll2004h . and are associated 
with AGN activity. 



3.2 Source counterparts at other wavelengths 

The NASA Extragalactic Databas^ (NED) was searched to find 
existing information on the IFRSs. The eight sources with the great- 
est 1.4-GHz flux density were present in the 1.4-GHz FIRST sur- 
vey. Four sourc es were present in t he 1.4-GHz NRAO VLA Sky 
Survey fNVSS: ICondon et alJll998h . a nd four were present i n the 
1.4-G Hz survey of the xFLS region bv lMorganti et alj j2004h . We 
use the lCondon et al.l J2003h 1.4-GHz data in preference to any of 
these surveys, due to the combination of high sensitivity and res- 
olution. One source (FLSGMRT J171436.6+594456) has a radio 
counterpart wi thin the 325-MHz We sterbork Northern Sky Sur- 
vey (WENSS: iRengelink et alj|l997l) . The flux density given in 
NED is 27.0 ±4.9 mjy, which is in agreement with the value of 



^ http://nedwww.ipac.caltech.edu/ 



30.8 ±0.4 mJy obtained by extrapolating the 610-MHz flux den- 
sity and spectral index (see section [33t . 

Eight of the sources had an optical coun terpart within 3 arcsec 
in the iJ-band images of lFaddaetij]j2004h . The i?-band magni- 
tudes (Vega) are listed in Table [T] where appropriate, and are close 
to the 50 per cent completeness limit of the observations, which 
was estimated tobe R = 24.5. Where no detection was found, this 
completeness level was taken as the upper limit on /?-band flux. 
No other optical or infrared counterparts were found for any of the 
IFRSs. 



3.3 Radio spectral index 

One of the few tools available for determining the type of IFRS 
being observed is their radio spectral index, due to the lack of de- 
tections in other wavebands. The spectral index of each source is 
listed in Table [T] and appear to be uniformly distributed between 
and 1.4 - a one-sample KS test gives a 99 per cent chance that 
the data are consistent with a uniform distribution, although the low 
number of sources limits the significance of this conclusion. 

We can use the spectral index as a discriminant to determine 
between different potential source types which make up the IFRS 
population. The lack of infrared detections, and VLBI detection 
of a compact radio core in one source indicate that obscured star- 
forming galaxies are unlikely to be the dominant source population 
- this is reinforced by the fact that only six of the sources have a 
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Figure 5. The luminosity - linear size plot for sources in the 3CRR cata- 
logue, along with the tracks that each of the IFRSs would make across the 
plot if they were located at redshifts between z = 1 (lower right) and z = 8 
(upper left). Each line represents the track of one IFRS, and the linear sizes 
given by the IFRS tracks are upper limits. Sources have been separated into 
flat-spectrum (a < 0.5; diagonal crosses for the 3CRR sources, solid lines 
for the IFRSs), steep-spectrum (0.5 < a < 1; dots for the 3CRR sources, 
dashed lines for the IFRSs) and USS sources (a > 1; upright crosses for the 
3CRR sources, dotted lines for the IFRSs). The flat-spectrum source 3C273 
is marked by the solid square. The horizontal line marks the FRI/FRII divide 
at 1.4 GHz. 



spectral index in the range 0.5 < a < 1, close t o the typical v alue 
for 'normal' star-forming galaxies of --^ 0.8 (e.g. ICondonlll992h . 

Three of the IFRSs are Ultra Steep Spectrum (USS) sources, 
with a > 1. Selecting sources with ultra-steep spectral in- 
dices is one method for efficiently finding high-redshift radio 
galaxies (HzRGS, e . g. Tielens et al. 1979; Rottgering et al. 19941 ; 



iBlundell et all Il998l ; Ide Breuck et al.l |2000|; iKlamer et alj 12006) . 
There are a number of mechanisms proposed to explain the re- 
lationship seen between spect r al ind ex and redshift (the z-a cor- 
relation) - see iKlamer et al. I 12006") for a recent discussion of 
these - with the high-redshift sources selected using this technique 
typically being compa ct, luminous Fanaroff-Riley Type II (FRII; 
lFanaroff&Rilevlll974h radio galaxies. 

The remaining five sources have flat spectra, with a < 0.5. A 
flat-spectrum source with compact morphology (see section|4Tl( is 
an indication of synchrotron self absorption occurring either due 
to core-dominated sources, or sources with a number of overlap- 
ping optically-thick reg ions aligned close to our line of sight (e.g. 
Ijarvis & McLur3l2002i) . 



4 MODELLING THE IFRS POPULATION 
4.1 Linear size constraints 

We list the deconvolved angular diameters of e a ch so urce in Ta- 
ble [T] taken from the catalogue of ICondon et alj ( |2003[) . From the 
larger diameter, we calculate upper limits on the linear size of the 
IFRSs at a given redshift, using the relationship between angular 
diameter-distance and z- We model the luminosity and maximum 



linear size for each source at redshifts between and 8, using the 
observed radio flux densities, and fc-correcting each source to ob- 
tain a rest-frame 1.4-GHz luminosity, Li 4, at each redshift, with a 
/^-correction factor of (1 (see e.g. Paper II). Fig.[5]shows the 

tracks of maximum linear size against luminosity for the IFRSs on 
a luminosity-size (P-D) diagram, for 1 < z < 8, along with the loca- 
tions of sour ces in the Revised Revised Thi rd Cambridge (3CRR) 
catalogue o f iLaing, Riley & Longaid Jl98^R the definitive list of 
bright radio sources in the northern sky at 178 MHz. The 3CRR 
catalogue has been chosen for comparison to the IFRS population 
due to the large amounts of photometry which are available on a 
significant number of bright radio sources. The break luminosity 
between Fanaroff-Riley Class I and Class II sources (FRI/FRII) is 
shown, adjusted to a 1.4-GHz value of 3.84 x lO^'* W Hz"' sr"' 
using a = 0.8. The flat, steep and USS sources are shown sepa- 
rately. 

The typical upper limit on linear size for the IFRSs is ~ 20 kpc 
at any redshift, and it is clear from Fig.[5]that the majority of 3CRR 
sources are too large to represent local versions of the IFRSs. In 
the following sections we model each class of source in turn, us- 
ing spectral energy distributions of suitable sources from the 3CRR 
catalogue in order to demonstrate that IFRSs can be modelled as be- 
ing high-redshift versions of well-studied local sources. Through- 
out this modelling, we fix the luminosity of the template 3CRR 
sources to match the radio flux density constraints for each IFRS, 
and use sources from Table[T]which have not been classified as hav- 
ing possible infrared counterparts where possible, since they are the 
more extreme examples of IFRSs. 



4.2 Modelling the flat-spectrum sources 

iNorris et al. I j2007l) claimed that, based on their 1.4-GHz detec- 
tion and 3.6-^m limit, the IFRS population may be represented 
by a high-redshift version of 3 C273. This source is a radio-loud 
quasar located at z = 0.158339 IStrauss et al j|l992|), with a single 
radio jet visible with linear s ize ~ 30 kpc H ahcall et al. 199 5|), L4- 
GHz flux density of 45 Jy jKellermann, Pauliny-Toth & Williamsl 
1 19691) . and spectral index of ~ 0. 1, giving it a luminosity of 
2.1 X 10'^ W Hz"' sr"'. It is not actually in the 3CRR cata- 
logue due to its low declination, but has similar properties to the 
flat-spectrum IFRSs (although it is slightly larger than the linear 
size constraints from Fig.[5]l- We have a numb er of additional con - 
straints on the spectra of IFRSs compared with lNorris et alj ( |2007|) : 



(i) We have radio detections of 1171453.8+594329 (the IFRS 
which fitted the spectra of 3C273 most successfully) at both 
610 MHz and 1.4 GHz. The luminosity of 3C273 was reduced by a 
constant (redshift-dependent) factor in order to match these obser- 
vations. 

(ii) We required that the 1 .4-GHz luminosity of the reduced ver- 
sion of 3C273 was great enough to remain above the FRI/FRII 
break luminosity. If the luminosity was to be reduced below this 
level, the FRII SED of 3C273 would be modelling a radio galaxy 
with the luminosity of an FRI source, which would be inappropri- 
ate. 

(iii) We place infrared upper limits on the IFRS population in the 
IRAC and MIPS bands, between observation-frame wavelengths of 
3.6 and 160 /im. 

(iv) We have optical detections, or upper limits, in the /?-band, 
which further constrains the redshift of each source. 

' http://3crr.extragalactic.info/ 
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Figure 6. The spectral energy distribution of 3C273, as it is observed (black 
line), and as it would appear if 3C273 was placed at z = 5 and reduced in 
luminosity by a factor of 15 (grey dashed line). The detections and upper 
limits on flux density for the IFRS J171453. 8+594329 are shown (diagonal 
crosses and upper limits) - the SED of 3C273 would be below the infrared 
detection limits, and match the optical photometry when placed at z = 5. 
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Figure 7. The spectral energy distribution of 3C305, as it is observed (black 
line), and as it would appear if 3C305 was placed at z = 1 and reduced in 
luminosity by a factor of 1.5 (grey dashed line). The detections and upper 
hmits on flux density for the IFRS J172315. 5+590919 are shown (diago- 
nal crosses and upper limits) - the FRI SED of 3C305 is unable to fit the 
infrared and optical data successfully. 



We obtained photometry on 3C273 from NED for compari- 
son with the IFRS data, and placed it at redshifts between and 
8, to examine how the SED that would be observed changes with 
redshift. Fig. |6] shows the observed SED of 3C273, and the SED 
that would be observed if 3C273 was placed at z = 5, and re- 
duced in luminosity by a factor of 15 (to a 1.4-GHz luminosity 
of 1.4 X 10^^ W Hz"' sr"'). This is a lower hmit to the redshift 
of the source (based upon fitting the optical/infrared photometry). 
Once the redshift is great enough that the SED shape can be fit- 
ted successfully, there is a degeneracy between the redshift and 
luminosity of the source, with the luminosity of 3C273 requiring 
a smaller amount of scaling at greater redshifts. The other flat- 
spectrum IFRSs can also be modelled with the SED of 3C273, with 
similar lower limits being placed on their redshift. The compact 
flat-spectrum sources in the 3CRR catalogue were also tested to 
see whether they could successfully model the flat-spectrum IFRS 
population, but less photometry was available on these sources, and 
none of them was as good a fit to the data as 3C273. Since the 
3CRR catalogue contains the brightest (low-frequency) sources in 
the sky, it should not be suprising that the flat-spectrum sources 
contained within it need to have their luminosity reduced signifi- 
cantly in order to fit the IFRSs. In carrying out the modelling in this 
section, we have assumed that a reduction in luminosity by a fac- 
tor of 15 does not fundamentally alter the properties of the quasar 
being modelled - it would be preferable to use lower-luminosity 
sources for this modelling, however there are few radio sources in 
the sky that have been as well-studied at optical and infrared wave- 
lengths as the 3CRR catalogue. 

4.3 Modelling the steep-spectrum sources 

There are six steep-spectrum sources with 0.5 < a < 1 within the 
IFRS sample, which follow very similar tracks in the P-D diagram. 
The 3CRR catalogue contains several local sources with steep spec- 
tra and linear size <20 kpc, with both FRI and FRII sources satis- 
fying the linear size requirements, and it is necessary to test both 
forms of SED in order to determine which class of source best rep- 
resents the IFRS population. 

Fig. [7] sh ows the spectrum of 3C305, an FRI radio galaxy at 
z = 0.041639 dMilleret al.ll2002l) . which has been placed at z = 1, 
and reduced in luminosity slightly (by a factor of 1 .5). While the ra- 
dio detections and far-infrared flux limits can be matched, the mid- 
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Figure 8. The spectral energy distribution of 3C67, as it is observed (black 
line), and as it would appeal' if 3C67 was placed at z = 4 and reduced in 
luminosity by a factor of 4.2 (grey dashed line). The detections and upper 
limits on flux density for the IFRS J172315. 5+590919 are shown (diagonal 
crosses and upper limits) - the FRII SED of 3C67 is able to fit the infrared 
and optical data successfully. 



infrared and optical hmits placed on the IFRS J172315.5+590919 
are inconsistent with the flux that would be observed from this 
source. No improvement was found to be possible at other red- 
shifts. In contrast to this. Fig. [U shows the spectrum of 3C67, an 
FRII radio galaxy at z = 0.3102 l lHewitt & Burbidge 1991), which 
has been placed at z = 4 and reduced in luminosity by a factor of 
4.2. All of the available photometry can now be fitted successfully, 
demonstrating that the steep-spectrum sources are better fitted by 
FRII templates than by FRI templates. A successful fit is also pos- 
sible with 3C67 placed at lower redshifts (i.e. 3C67 at z = 3 would 
require a reduction in the total luminosity by a factor of ~ 9 to fit 
the radio data), but the origin of the optical and radio luminosities 
in FRII galaxies is very different and this amount of luminosity 
scaling may not be justified. The lower limit on redshift for this 
source is z = 2.5, placed by requiring the 1.4-GHz radio luminosity 
to remain above the FRI/FRII break. 

4.4 Modelling the ultra-steep spectrum sources 

There are three sources within the IFRS sample with a > 1. One 
of these, J171928.4+583927, has a spectral index of 1.04±0.09 
and is better fitted by the SED of 3C67 than to any of the 3CRR 
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Figure 9. The spectral energy distribution of 3C186, as it is observed (black 
line), and as it would appear if 3C 1 86 was placed at z = 2 and reduced in 
luminosity by a factor of 90 (grey dashed line). The detections and upper 
limits on flux density for the IFRS J171315. 5+590302 are shown (diagonal 
crosses and limits) - the SED of 3C186 is able to fit the radio and opti- 
cal data, and the infrared limits at 3.6 and 4.5 ilm (the only points where 
photometry on 3C186 is available). 



USS sources, while the other two have spectral indices of 1.37 
and 1.38. The 3CRR catalogue only contains two compact sources 
with ultra- steep spectra - of the two, 3C186 (an FRII source at 
z = 1.063: lLvndsetai]|l966h has more available photometry, and 
Fig. |9] demonstrates that it is possible to fit the radio and optical 
data for IFRS J171315.5+590302, along with the 3.6- /im IRAC in- 
frared limit, by placing 3C 1 86 at z = 2 and reducing its luminosity 
by a factor of 90. As seen in previous sections, placing the source at 
higher redshift allows this luminosity scaling to be reduced while 
still being able to fit the observed photometry - at z = 4 the lumi- 
nosity of 3C186 needs to be reduced by a factor of ~ 15, while at 
z = 6 the luminosity only needs to be reduced by a factor of --^ 5. 



4.5 Comparison with the 'Optically Invisible' radio sources 
of Higdon et al. (2005) 

iHigdon et alj ( |2005[) describe a population of sources which are 
radio-bright but 'optically invisible', with the majority of their 
optically-invisible radio sources (OlRSs) also being undetected in 
the mid-infrared with MIPS. While the majority of their OIRSs are 
not particularly radio-bright, there are two sources with 1.4-GHz 
flux density > 1 mjy, with no detections in optical images having 
Bw ~ 27.0, R ~ 25.7, 7 ~ 25.0 (Vee a magn itudes), or at 24 ^um, 
with 524 < 0.3 mj}0. iHigdon et aU \2Q0§ ) conclude that these 
sources are powered by AGN activity (principally based upon the 
non-detect ions at 24 jxm ) . IRA C information is added to the OIRS 
sample by IHigdon et alj fcoOSi) . who find that at least 34 per cent 
of the sources do not have detections at 3.6 /im. They conclude 
that the OIRSs which are undetected by IRAC appear to represent 
a population of powerful radio galaxies, at z > 2. 

The Infrared-Faint and Optically-Invisible radio source popu- 
lations appear to be similar in nature, although with slightly differ- 
ent selection effects. While it remains to be confirmed that these 
sources are high-redshift radio galaxies (potentially through ex- 
tremely deep multi-wavelength optical and infrared observations). 



There are two further radio sources with S1.4 > 1 mJy and no optical 
detections, outside the coverage area of MIPS, and one further source which 
is desc ribed as an O IRS in Hi gdon et al. 1 .2005.) . but no longer thought to be 
one bv lHigdonet"alJ<2008i) . 



it is not necessary to require new, exotic source types in order to 
explain the IFRS or OIRS population. 



5 CONCLUSIONS 

We confirm the existence of a class of bright radio sources which 
are undetected in deep infrared images. These sources are rare - 
we find a total of 14 sources with 5i.4 > 0.5 mJy and no clear 
infrared counterparts. Ten of these have a flux density of greater 
than 1 mJy, a source density of roughly 2.5 deg^^, and compa- 
rable t o th e source density of brigh t IFRSs found by iNorris et al.l 
( I2OO6I) and lMiddelberg et al.l |20o9). There are eight sources with 
5i.4 > 0.5 mJy that show no evidence for possible infrared coun- 
terparts in the deep IRAC images of the Spitzer extragalactic First 
Look Survey field. 

We use source stacking to place upper limits on the mid- 
infrared flux densities of these sources, and use either detections 
or upper limits in the i?-band to further constrain the SED of 
each source. By using multi-frequency radio data, we demonstrate 
that the sources vary in type between flat-spectrum quasars with 
a --^ 0.1, and USS radio sources with a ~ 1.4, and should not be 
thought of as a single source population. 

The possibility that the Infrared-Faint Radio Source popula- 
tion is made up of high-redshift luminous radio galaxies appears 
increasingly likely - the lack of detection in infrared images would 
seem to rule out highly obscured star-forming galaxies, while the 
existence of flat-spectrum sources implies that at least some of the 
sources have an AGN origin. We place upper limits on their over- 
all linear size through the k nowledge of their m aximum decon- 
volved angular size from the ICondon et alj ( |2003[) 1.4-GHz cata- 
logue, and show that they must be < 20 kpc at any redshift, making 
them smaller than the majority of bright radio sources in the local 
Universe, as given by the 3CRR catalogue. By looking at the lo- 
cations of these sources in a P-D diagram, for a range of redshifts, 
we identify bright local radio sources with similar characteristics, 
and demonstrate that all three classes of IFRS can be modelled as 
less-luminous versions of well-studied Fanaroff-Riley Type II radio 
galaxies placed at high redshift. 
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